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Abstract. For the configuration of a sphere in front of a plane we calculate the first 
r~~"i, two terms of the asymptotic expansion for small separation of the Casimir force. We 

^3 ' consider both Dirichlet and Neumann boundary conditions. 

1. Introduction and discussion 

> ' 

^ ■ High precision measurements of the Casimir force are of remarkable actual interest. For 
instance, these measurements provide stronger constraints on hypothetical interactions 
in the sub-micrometer region pQ and give hope to measure thermal forces. The high 

CN 1 precision is achieved in a setup of measuring the force between a flat plate and a sphere. 

qq ■ The typical radius of the sphere is R ~ 200/xm (although recently [2] a sphere with 
R ~ 80/xm was used) and the separation ranges from contact till L m lfim. This setup 
provides a large area of opposing surfaces avoiding difficulties from keeping the surfaces 

^ ■ in parallel. The ratio of these two sizes, 

is a small parameter describing the deviation of the geometry from the plane parallel 
one. 

The calculation of the Casimir force for arbitrarily shaped surfaces was a long 
standing problem and one was left with unsatisfactory approximate methods. The most 
important of them is the proximity force approximation (PFA) known from [3]. The 
method takes the force density known from the plane parallel case and integrates it 
over the curved surface. Clearly, this methods works only for small deviations and 
from the very setup it was impossible to get higher order corrections or information on 
the precision of the approximation. Nevertheless, PFA worked remarkably well in the 
comparison of the precision measurements [U [2] with theoretical expectations. Along 
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with this, from increasing the precision of the measurements there appears a call to go 
beyond [2]. 

The evaluation of the Casimir force is in general plagued by the ultraviolet 
divergences even what concerns the force between distinct bodies although these are 
finite. The point is that divergences are present in intermediate steps. One has 
to introduce some regularization, say in a mode sum representation, to subtract the 
divergent terms (even if these do not depend on the separation), and to perform the limit 
of removing the regularization. This procedure works well in simple geometries with 
separating variables in the corresponding wave equation. However, for a more general 
geometry, where one would have to calculate eigenvalues numerically, this hampers the 
calculations quite much. 

There exists approximate methods which gave important results. For instance, the 
semiclassical approximation in [5l E] has put the PFA on a solid base within quantum 
field theory. Also the optical path approximation [7] was very important in this context, 
since it also confirmed the leading-order PFA from a field-theory perspective and gave 
somewhat better results at next-to-leading order. An exceptional role play the so called 
'world line methods' [8], where a functional integral representation of the Greens function 
is computed numerically. Regrettably, this method is restricted to Dirichlet boundary 
conditions. 

Progress happened recently with the representation of the Casimir interaction 
energy between two bodies in terms of a functional determinant which is free of 
ultraviolet divergences [9], [10] . For any fixed distance between the bodies, the energy is 
represented by multiple convergent sums and integrals. These allow for direct numerical 
evaluation [9]. However, for small e like the value in ([1]), the computations become too 
complex and are restricted to e > 0.1. 

A solution of this problem was found in [11] where the first correction beyond PFA 
was calculated for a cylinder in front of a plane. Technically, an asymptotic expansion 
for small e of the representation in terms of a functional determinant was derived. In 
the leading order, the PFA was re-obtained and in the next-to-leading order the first 
correction beyond PFA was calculated analytically. This was done for Dirichlet and for 
Neumann boundary conditions, 
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The sum of these two, by virtue of the wave guide geometry, delivers at once the result 
for the electromagnetic case. It is remarkable, that the result for the Dirichlet boundary 
conditions was with good precision confirmed by the world line methods in [12J. 

In the present paper we use the method developed in [TT] to calculate the Casimir 
interaction energy for a sphere in front of a plane. We perform the calculation for 
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Dirichlet and for Neumann boundary conditions. The result is 




(3) 



With these formulas we re-confirm the PFA in leading order. In next-to-leading order 
we obtained analytically the first correction beyond PFA. Again, as in the cylindrical 
case, the Dirichlet case is confirmed by the 'world line methods' calculation in [12J. 
Since in this geometry the polarizations of the electromagnetic field do not separate, ([3]) 
does not include the result for the electromagnetic case which is left as a future work. A 
general approach for that was considered in [T3] along with specific calculations of the 
Casimir interaction between spheres over a wide range of separations, including rather 
short distances. However, if expecting the corresponding result to be of the same order 
as ([3]) one can conclude that these corrections do not affect the current high precision 
measurements of the Casimir force but must be taken into account in future ones. 

In the remaining part of this paper we follow [TT] and represent the calculations 
resulting in ([3]). The essentially new element is a certain asymptotic expansion of the 
Clebsch-Gordan coefficients which was not found in literature. 

2. Asymptotic expansion of the Casimir energy for small separation 

The expression of the Casimir energy for a sphere in front of a plane in terms of a 
convergent functional determinant is known from [101 E] ■ Here we use the notations 
adopted in [TT]. The energy under consideration is given by 



Here the trace is over the orbital momenta I and m in a representation in orbital 
momentum basis of the matrix elements 



where I v and K v are modified Bessel functions, a = L + R is the distance from the 
plane to the center of the sphere and the symbols in brackets are the 3j-symbols. It is 
convenient for the following to expand the logarithm in (j4j), 




(4) 
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For fixed R and a, the integral and all sums in this formula are convergent. With 
e = L/R, ([1]), we note a = R(l + e) and for small e the convergence gets lost. As 
for the frequency u this can be seen from the arguments of the Bessel functions. Using 
their asymptotic behavior for large argument, I v {z) ~ exp(z) and K u (z) ~ exp(-z), the 
integrand in the matrix element behaves like exp(— 2ue) and ceases to decrease for e = 0. 
Similar arguments hold for the sums involved. As a consequence, a simple expansion 
for small e does not work. 

Next we rewrite Eq.([6]) by making the substitution u> — > uj/R and by changing the 
summation indices in the products of the matrix elements, 



oo oo oo oo 



B =-^ihl ^ E EE -E 

s=0 JU m=-oo l=\m\h=-l l 3 =-l 

■ Ai,i+h^i+h,i+h ■ ■ ■ Ai + i Si i . (7) 

The strategy of the calculation for small e is as follows. We assume that the 
dominating contribution comes from the large values of the frequency u and all 
summation indexes involved in (JTj). We are interested in the asymptotic expansion 
for small e, hence we substitute all sums by integrations. The error introduced in this 
way is supposed to be exponentially small like in the difference between sum and integral 
in the Abel-Plana formula. After that we make the substitution of variables in (H), 

h = ^ ^m, . . . ,l 8 = \j-^ n s , (8) 

manifesting that the main contributions come from u ~ 1/e, from the 'main diagonal' 
matrix index / ~ 1/e, from the 'off-diagonal directions' U ~ {i = l,...,s) and 

from the magnetic quantum number m ~ 1/ \ft. The variable r is the cosine of the angle 
in the (u, /)-plane. After this we expand the matrix elements for e — >• and calculate 
the remaining integrals. 

In the new variables the expression for the energy reads 

e- 2 1 f°° , f l drr f°° da 
E ~ > dtt 1 ' 

Att R ^ c 4-1 /„ 

oo 
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M=(—\ \ l+h ■ . . . ■ A l+lsjl (10) 



where 



\ e 

collects the information from the matrix elements. 

In the appendix A we calculate the asymptotic expansion of the matrix elements (jSJ) 
for e — >• with the variables substituted according to Eq. (jSJ). Using Eq. fl38|) obtained 



Casimir force beyond PFA 



5 



there we get for the asymptotic expansion of the energy the expression 



e -(s+l)M 2 



^ ; ^...r^e-JM. (11! 



with 

s 

77i = ^ (rii - n i+1 ) 2 . (12) 

Here, we have formally to put uq = n s+ \ = 0. The non- leading information is contained 
in 



m™ = 1 + *fi E + £ E «S3„<£» + E <!».+, + ■ ■ 



i=0 \0<i<j'<s i=0 

where the sums result from the products of the matrix elements in (11 01) . Now it is 
possible to carry out the integrations. This can be done in complete analogy to the 
cylindrical case. For instance, in the leading order we note 

dm f°° dn„ 1 . , 

' e~ Vl = - . (13) 



^ " ' J-00 y/s + l 

Another square root of (s + 1) comes from the integration over \x and two more inverse 
powers of (s + 1) come from the integration over t. Collecting all factors together 
we immediately come to the expression know from the PFA thus reconfirming that 
approximation independently. For the next-to-leading orders we also use the same 
techniques as in the cylindrical case. Thus the integrations over the rij are Gaussian 
integrations with a somehow involved combinatorics. These are completely described in 
the appendix C in [lljjl After that the integrations over the remaining variables can be 
carried out too. All these integrations are either Gaussian or simple exponentials with 
polynomial factors in front. However, the number of terms in the intermediate steps is 
quite large. In fact, the calculations were performed machined using to a large extend 
the same scripts which served in the cylindrical case. The result is the remarkably simple 
coefficient 1/3 in Eq.(j3]). 

Finally, we describe the changes which come in if using Neumann boundary 
conditions instead of Dirichlet ones. First of all, from the Neumann condition on the 
plane we have a sign change in the expression for the energy which now reads 
1 f°° 

E=— dcu tr,, m In (<?,,,, + A%,) . (14) 



The derivatives on the sphere result in changed matrix elements Af v . From the 
derivation of the matrix elements it follows that the radial arguments of the matrix 

% Regrettably, that appendix contains a number of misprints. These are corrected in the electronic 
preprint, |arXiv:hep-th/060 2295v2. 
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elements r and r' are just in the arguments of the Bessel functions I v+ \(ujr') and 
in ([3]) before putting them on the sphere, i.e., before putting r = R and 
The derivative in the Neumann boundary conditions acts just on these 



r' = R. 



arguments. Consequently, we have to substitute I v+ i(uor') — > (rj^ + u) Iy + i(ur') where 
u is the parameter which appears in Robin boundary conditions. In fact, it is possible 
to formulate different types of Neumann conditions by multiplying the wave function 
with some power of the radius, which by means of r<9 r r u $ = r u (rd r + m)$ is equivalent 
to a Robin condition. For a single sphere, one has to put u = 1/2 in the electromagnetic 
case. 

Now we discuss the changes which happen in the asymptotic expansions of these two 
Bessel function in the matrix element. It can be easily verified that for the asymptotic 
expansion of the Bessel functions the following formulas hold, 

I u = (rd r + u) I v {yz) = vzl' v {uz) + uI v (pz), (15) 

K u = (rd r + u) K v (vz) = vzK' v (vz) + uK u (uz). 

Using the asymptotic expansions these expressions can be rewritten in the form 



/ 



Wl + z 2 Uuz) 



u 



+ 



fc>0 

"[ft 

k>0 



V 



Vl + z 2 K u (uz) 



u 



V 



+ 



/ 

E(-i) fc ^\ 

fc>0 



\k u k 



(16) 



where Uk and Vk are the known Debye polynomials. Now, to leading order for large u, 
we get for the expressions in the brackets 

u 



zVl + z 2 I v {pz) 



Vx — Ui 

1 + — — - + 



V 



Vl + z 2 K v 



vz) 



uVT+z 2 

Ui u 



17) 



V V\f\ + z 

These expressions have to be inserted into (jSJ). It is clear that the factors vyl + z 2 
cancel and that the expansion in the brackets becomes 



1 + 



+ 



uVl+z 



1 



°l-s 



The dependence on u drops out and the remaining factor is the correction factor which 
must be inserted into the calculation of the Dirichlet case. We note that the sign 
coming from the Bessel function K v [yz) in ( [16]) compensates the changed sign in ( TH1) . 
These are the only changes which come in and performing the corresponding calculation 
one comes to the second line in ([3]). In this way the first correction beyond PFA for 
Neumann boundary conditions is calculated. From the above discussion it is evident 
that in leading order the same result as for Dirichlet boundary conditions appears thus 
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re-confirming again the PFA which is insensitive to the boundary condition and that 
the first correction beyond does not depend on the type of Neumann condition, i.e., it 
does not depend on the parameter u. 

This work was supported by the research funding from the ECs Sixth Framework 
Programme within the STRP project "PARNASS" (NMP4-CT-2005-01707). 
V.N. was supported by the Swedish Research Council (Vetenskapsradet), grant 621- 
2006-3046. 

Appendix 

In this appendix we calculate the asymptotic expansion of the matrix elements Ai t i>, (JSJ). 
We start with the substitution I" = I + I' — 1v of the summation index, 



substituted according to (jHJ). To obtain the asymptotic expansion of Ai + i it i + i i+1 we insert 
the asymptotic expansions for the Bessel functions and for the 3j-symbols and extend the 
sum over v till infinity. We face with the necessity to derive the asymptotic expansion 
of the 3j-symbols, whereas the asymptotic expansion for the Bessel functions is well 
known. 

In general, asymptotic expansions of the 3j-symbols, or of the related by means of 



Clebsch-Gordan coefficients, are well investigated. However, the main attentions is paid 
usually to semiclassical expansions aiming for computational tools for large quantum 
numbers. A comprehensive and modern representation is given in [13]. Unfortunately, 
it does not cover our case. The point is that [H] considers the limit where all 
quantum numbers become large whereas in our case the magnetic quantum number 
m grows slower that the other ones. Nevertheless [13] gives a very convenient integral 
representation which serves as the starting point in our case too. We used Eq.(2.11) in 
[14] which in our notations reads 




We are interested in the asymptotic expansion of Ai + i u i + i i+1 for e — > with the variables 




(20) 




c; 



d",0 

l,m;l' ,—m 



(_4)P"+W)/2jVj. 



l,m;l'-m;l",0 




■ (sm(6 - 0)) 



l"+l+l' 



d9d(j) 



(21) 
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whereby we shifted the angles 9 and by ir/2 which resulted in the factors cos# and 
cos0 in place to the factors sin9 and sin</> in [T3] . The factor iVj m; ^ - m -i» ,q is defined 
by 

Ni im ;l',-m;l",0 = V + 2n^ + 2n m - 4z/ + 1) 

(/ + rii + m)\{l + rii — m)\{l + n i+1 + m)!(Z + — rn)\ 
(41 + 2m + 2n i+l -2v + l)!(2z/)! 
(2/ + ni + n i+1 -2z/)! 



((2/ + 2rii - 2z/)!(2Z + 2n i+1 - 2z/)!) 3 

in our notations. 

The asymptotic region we are interested in is determined by Eq. ([8]) . So we have to 
use I" = l + V — 2z/, I — * l + li and /' — > l + k + i in Eqs. fl2~TT) and fl22|) to perform calculations 
in the asymptotic region. First of all we consider iV/ jm; // _ m; z// j0 - Its asymptotic expansion 
follows simply from Stirlings formula. For small e it holds 

\j2ixa 



(22) 



a b 
- + —= + c 
e sje 



cxp 



'.2 1 



-1 + In - ) + — In - + cln - + 

Ve e e 



2a 



1 + 



b (a(6c + 3) - b 2 ) y/l 
6a 2 



+ 



(o 6 - 12aco 4 + 9a 2 (4c 2 - 1) b 2 + 6a 3 (6c 2 + 6c + 1)) e 



72a 4 



+ O (e 3 / 2 ) 

Now we make the substitutions ([8]) in Ni )m; i> - m -i" ,o? (|22l) . and obtain with (|23l) 



^ r Z,m;/',-m;i",0 

where the function 



7T 3 / 4 



2-H-2{li+U+i)+2u+l/A y+3/4 e m 2 /l jy/^ 



(23) 



(24) 



ft, \ i ((n, + n w )(4^ 2 -4z/-3)) 

iV(e) = 1 4(^5 ^ 

+96(1 - 4/i 2 + 4i/)nin i+ i 

+ (16/i 4 - 48/i 2 - 3 {Su 2 + Av- 5)) r) e 

+ O (e 3 / 2 ) (25) 

collects the subleading contributions. 

Next we calculate the integral in fl2T|) using the saddle point method. The 
considered integral I we rewrite in the form 

r-Tr/2 ptt/2 



\2v 2im(e-<j>) 



-{l"+l-l')h{<f>)-(l"-l+l' 



tt/2 J -it/2 



with 



fc(0) =lncos0, g( 



sin w — 



2u 



(26) 



(27) 
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For small e, because of I" + I — V = 2tr/t + \jAt/t{n,i — n i+ i) — 2v, the dominating 
contributions come from the minimum of the functions h which are reached in 9 = 
and = 0. Therefore we expand these functions using h(x) = x 2 /2 + x 4 /12 + . . .. The 
leading order contribution in the exponential is then exp(-^(fl 2 + 2 )). Taking into 
account that the function g depends on the difference, 9 — 0, we are motivated to make 
the substitution of variables 



9 =Vs? K+ " ) ' W^-"'- (28) 

After that the integral can be rewritten now in the form 

(r> \ 1+f />oo poo 
A / / y^e^^m^e-^d^drj, (29) 
J J — OD J — OO 

In the sense of the asymptotic expansion we are interested in we stretched the integration 
until infinity. In (1291) we defined the function 

g(£, v ) = g(9, 0) e Hi"+i-i'M4>)-(i''-i+i'M0)+e+v\ (30 ) 

Obviously, it has an expansion in powers of y/e, 

<?(£> V) = 1 + (-(n< + n i+1 )r] 2 + 2(n i+1 - m)£r] - (m + n i+1 )£ 2 ) 

y/tr 

+ ((6(n< + n i+1 f -t) V 4 + 24 (n 2 - n i+l 2 ) £ V 3 
+ (6 (Qn 2 - Am+im + Qn i+1 2 - r) £ 2 + Apt) rj 2 
+24 (n 2 - n l+1 2 ) fr] 

+ (6(n, + n m ) 2 - r) e 4 + 12^ 2 ) ^ + O (e 3 / 2 ) . (31) 

Now we carry out the integration over £. It is Gaussian and delivers a new function 
<?i(e,7/), 

with the property pi(0, r/) = 1. Finally, we take the expansion (1241) of the prefactor 
and the expansion ([3"2"j) of the integral together and insert into Eq. f[2~Tj) for the Clebsch- 
Gordan coefficients. This gives 

cK,-,- m = ^ 2-«/4VM rl /4 ft(£i „) t (33) 

where the new function h(e, //) collects the sub-leading contributions, 

o2i/..| /"oo 

h(e, n) = ==-. — / r] 2u g 2 {e, 77) f&^^f+f dq . (34) 



rt 1 (e,r ? )e- ,?+W2T ^dr/ (32) 



—00 



v^(2^)!i-c 

Here we defined g 2 (e,ri) = ni(e)gi(e,T]) which must be re-expanded in powers of \/e 
too. The function h(e,fi) has the property h(0,fi) = 1. It is possible to obtain a more 
explicit representation of this function by carrying out the integration. However, for the 
application below this proved to be not useful. 
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The calculation of the asymptotic expansion (1551) of the Clebsch-Gordan coefficients 
was the most difficult task. It remains to insert them into the matrix elements Ai i>, 
ffT9]) . These symbols enter twice. The second time they go with zero magnetic quantum 
number m. It is possible to put m = directly in the asymptotic formula ( 133]) . In ( |34l ) 
then the integration can be carried out explicitly. Putting these formulas together with 
the well known asymptotic expansion of the Bessel functions (which is displayed, for 
instance, in Eq.(B.l) in [11]) and we come to 

re'"" r dr, 2+WI „ + „> ^ rf" /1-tV 

where C g collects the sub-leading factors. Its dependence on v is polynomial. Hence, 
the sum over v can be carried out, 

^" C. = S&C„ (36) 



^ v\ V 1 + r 



where now C a collects the non-leading contributions. The complete expression is too 
lengthy to be shown here. So we restrict ourselves to show the order ^/e, 

C a = 1 + Q r ) A 2r\r + l)n\ - 2n l+1 r 2 (r + l)n 2 
2\jtT{j + 1) v 

+ {{-2n 2 i+1 - At + 1) r 3 + (V - 2n 2 +1 - 4t) r 2 - 8r/ 2 r 

+r - 2/2 2 (t + 1) + 2) n, + ( (2n 2 +1 - 4t - 3) r 3 (37) 

+2 (2?7 2 + n 2 +1 - 2t - 2) r 2 - 8r/ 2 r + r - 2/i 2 (r + 1) + 2) ) J~e. 



With (1361) the integration over rj in fl35|) becomes Gaussian and can be carried out. 
Together with the remaining factors the result collects into the final form of the matrix 
elements up to the e order, 



A l+k}l+k+1 = Va3 e ^ (l + a2(n;,n m )Ve + ^(rii,n i+ i)e) , (38) 

where 

r/ as = -2t + (rii - n i+1 ) 2 (39) 
is the same factor as in the cylindrical case [TT] and the non-leading coefficients are 



4 1 

2^r 2 



<,n i+1 = :r ^ 7 (2n 3 r 3 - 2n 2 n l+1 r 3 



a ni,n i+1 



+ m {{-2n 2 l+1 - At + 1) r 3 - 2/i 2 (r 2 - 2)) 

+ n i+1 ( (2n 2 +1 - At - 3) r 3 - 2^ 2 (r 2 - 2) ) ) (40) 

^ (4 ( 3 r 4 - 9r 2 + 6n 2 (r 2 - 2) 2 + 6n 2 +1 (r 2 - 2) 2 

+12^^+1 (r 2 - 2) 2 + 4) ^ - 12r (4r 2 (r 2 - 2) 
-2 (4 (n 2 +1 +t)r 4 + (-8n 2 +1 - St + 5) r 2 - 8) 



n 2 
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-4n i+ i ((4* + 2)r 4 + (1 - 8*)r 2 - 4) m - 4*r 4 - r 4 + 2r 3 
+4*r 2 - 3r 2 - 4r + 4n 4 +1 r 2 (r 2 - 2) 

(~8(t + l)^ 4 + 2(8* + 3)r 2 + 16) + 4) fi 2 
+ r 2 (3 (8nf - 16n m nf - 4 (2n 2 +1 + 8* + 5) n 4 
+ 16n m (2n 2 +1 - l) n 3 + (-8n 4 +1 + 8(8* + 5)n 2 +l 
+32* 2 - 22) n 2 

+4n m (-4n 4 +1 + 12n 2 +1 + 16t 2 + 16* + l) 
+8nf +1 + 16* 2 + 8* - 4n 4 +1 (8* + 13) 

+n 2 +1 (32* 2 + 64* + 58) - 5) r 4 + 12 (2n 2 - 4n m n; + 2n 2 +1 
-4* - 1) r 3 + (28n 4 - 16n m n 3 - 12 (2n 2 +1 + 4* - 3) n 2 
-8n i+ i (2n 2 +1 + 12* + 3) n» + 28n 4 +1 - 24* 

-12n 2 +1 (4* + 5) + 9) r 2 + 12) ) . (41) 
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